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Abstract-Per\sta\\\c pumping of a non-Newtinian 
micropolar fluid in an inclined channel is studied. 
The analysis has been carried out in the wave frame 
of reference with long wavelength and zero 
Reynolds number assumptions. The velocity 
distribution and stream functions are obtained. The 
pressure rise (drop) over one wavelength is 
obtained. The velocity field, the stream function, , 
the volume flux and the pressure rise over one cycle 
of the wave and frictional force are obtained and on 
the flow quantities are discussed graphically. 
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I. Introduction 

Peristaltic pumping is a form of fluid transport, 
generally from a region of lower to higher pressure, by 
means of a progressive wave of area contraction or 
expansion which propagates along the length of a tube- 
like structure. Peristalsis occurs naturally as a means of 
pumping biofluids from one place of the body to another. 
This mechanism occurs in the gastrointestinal, urinary 
and reproductive tracts and many other glandular ducts 
in the living body. The early reviews of Ramachandra 
Rao and Usha [1] , Jaffrin and Shapiro[2] , Manton [3], 
Brasseur et al. [4], Srivastava and Srivastava [5], 
Provost and Schwarz [6], Shukla and Gupta [7], Misra 
and Pandey [8] ,Rao and Rajagopal [9],Kavitha et.al[10] 
Vajravelu et al. [11-15], Subba Reddy[16&17], Srinivas 
[18,19] deal with the peristaltic transport of viscous 
fluids through tubes and channels having impermeable 
flexible walls. 



In classical continuum theory a body is assumed to 
be a dense collection of point masses in which there is 
no internal structure. In the motion of a volume element 

Av it is assumed that the individual motions of material 
points coincide with the motion of centre of mass of the 
volume element Av. In this case the density p of the 

volume element Av is independent of the size of Av 
and independent on it's location in space and the time t. 
Eringen [20] reported that this is not true 

asAv^O.The density p shows an increasing 

dependence on the size of theAv, when Av is less 

than a critical value Av* . Classical continuum theory 
cannot explain the mechanical behavior of rheologically 
complex fluids, such as liquid crystals, colloidal fluids 
and blood. Due to this fact a new approach was 
necessitated. There are several approaches to the 
formulation of microcontinuum theories of fluids such as 
simple deformable directed fluids, dipolar fluids, polar 
fluids, simple micro-fluids, micropolar fluids, etc. All 
these consider the existence of couple stresses and 
body couples. 

Eringen [20,21] reported the theory of micropolar 
fluids in which the fluid micro elements undergo 
rotations without stretching. Micropolar fluids are 
superior to the Navier-Stokes fluids and they can sustain 
stresses and body couples. Here the micro particles in 

the volume Av rotate with an angular velocity about 
the centre of gravity of the volume in an average 
sense and is described by the micro rotation vector 

£1 .The micropolar fluids can support stress and body 
couples and find their applications in a special case of 
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fluid in which micro rotational motions are important. 
Airman and Cakmak [22] discussed three basic viscous 
flows of micropolar fluids. They are Couette and 
Poiseuille flows between two parallel plates and the 
problem of a rotating fluid with a free surface. The 
results obtained are compared with the results of the 
classical fluid mechanics. Srinivasacharya et al.[23] 
made a study on the peristaltic pumping of a micro polar 
fluid in a tube. The gravitational effects are also 
important in peristaltic pumping. In view of this, we have 
considered the peristaltic pumping of a micropolar fluid 
in an inclined channel. This mathematical model may 
be useful to have a better understanding of the 
physiological systems such as blood vessels. The 
velocity field, the stream function, the volume flow rate 
and the pressure rise are obtained and results are 
discussed through graphs. 



y 




Figure 1. Physical model 



1 1. Mathematical Formulation and Solution 

Consider the peristaltic pumping of a micropolar fluid 
in an inclined channel of half-width 'a'. A longitudinal 
train of progressive sinusoidal waves takes place on the 
upper and lower walls of the channel. For simplicity we 
restrict our discussion to the half-width of the channel 
as shown in figure. (1) 

The wall deformation is given by 

H(X, t) = a + bsin — (X-ct) (1) 

Where b is the amplitude, X is the wavelength and c is 
the wave speed. 



Equations of Motion 

Under the assumption that the channel length is 
an integral multiple of the wavelength X and the 

pressure difference across the ends of the channel is a 
constant, the flow becomes steady in the wave frame 
(x, y) moving with velocity c away from the fixed 

(laboratory) frame (x, Y)- The transformation between 

these two frames is given by 

X= X-ct; y= Y; u(x, y)=U(X-ct,Y); v(x,y)=V(X-ct,Y) (2) 

Where U and V are velocity components in the 
laboratory frame and u, v are velocity components in the 
wave frame. In many physiological situations it is 
proved experimentally that the Reynolds number of the 
flow is very small. So, we assume that the flow is 
inertia-free. Further, we assume that the wavelength is 
infinite. 



Using the Non-dimensional Quantities. 

_ u _ x _ y _ pa 2 — (la , H 

u = — ; x = — ; y = — ; p = - ; £2 = ; h= — 

c A a kc[i c a 

The non-dimensional form of equations 
governing the motion (dropping the bars) is 
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d 2 u XT 3£2 /, XT x dp . 

+ N- (1 - N) ^+ r|sina = 0 



dy 2 dy 



2-N d 2 £l du 



dx 



m 



dy 2 dy 



-2£l = 0 



Where N = 



|Li + k 



is coupling number 



(3) 



(4) 



Where P = 



dp 

dx 



f 



£1 - A 2 coshmy + A 3 sinh my - 



V 



(l-N)p-r|sina 
(2-N) 



y-4 



J 



m 



using the boundary conditions (5) to (8) in 
(9) and (10),we obtain the velocity of the fluid and 
micro rotation velocity as 



£1 is the micro rotation velocity 

u is the velocity 

JLX is the viscosity of the fluid 

k is the micropolar viscosity 



m is the micropolar parameter 



u = 



(l-N)p-r|sina 



(2-N) 



N 2 
§ sunny — ooshny+y 

m 



-1(11) 



-2N 

where Sj = — — 



m 



Di 



p is the fluid pressure 



S 2 = 



2NDj 2ND 1 f sinhmh , Ncoshrrh % 2 



m S3 



V m 



my 



m 



The non-dimensional boundary conditions are 



du 

do. 

3y 



= 0 at y = 0 



= 0 at y = 0 



u = -1 



at y= h(x) 



£2 = 0 at y= h(x) 



Solution 



(5) 



(6) 

(7) 
(8) 



The general solution of (3) and (4) is given by 



u= — (Aj sinhny+Aj GOshny) + 
m 



^(l-N)P-r|sina^ 



V 



(2-N) 



y+A 



J 



m 



So = 2 cosh mh - N 



D, = mh - sinh mh 



(l-N)P-r|sina A 
£2=A 2 coshmy+A 3 sinhmy-- ^ y 



(2-N) 



m 

(12) 



where 



mN [( 1 - N) P - r) sin a] [ mh - sinh mh] 
1_ (2-N)(2coshmh-N) 



A. = 



2 [( 1 - N) P - r| sin a] [ mh - sinh mh] 



m(2-N)(2coshmh-N) 



(9) 



A 3 = 



(l-N)P-r|sina 



(2-N) 



m 
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Integrating the equation (11) and using the condition 

\|/ = 0 at y = 0 , we get the stream function as _ 

Q = q + 1 (16) 



(l-N)p-r|sina 



(2-N) 



SjCOshrry Nsinhny^y 3 
m m m 3 



(13) 



■y 



The Pumping Characteristics 

Integrating the equation (15) with respect to x 
over one wavelength, we get the pressure rise (drop) 
over one cycle of the wave as 



The volume flux q through each cross-section in 
the wave frame is given by 



n 

= judy 



(14) 



0 



4>=J 



1 



^cohnh Nsirhnh h 3 „ 
Sj + +S>h 

v m m m 3 j 



&+rpna 



where 



Si = 



-2N 



m 



q= 



(l-^P-r)sina 



(2-N) 



S,ooshrrh Nsinhrrh^h 
m m m 3 



-h 



S 2 = 



2NDj 2NDj ( sirfirrh h\ Ncoshnii l 2 



S3 



+ 



V m my 



m 



(17) 



The pressure rise required to produce zero 
average flow rate is denoted by AP 0 . Hence AP 0 is 



given by 



0 



1 



„ arfinti Nsihnh h ^ 

S| -^+S>h 

V m m m 3 j 



ck+rpna 



S a = 2 cosh mh - N 



D, = mh - sinh mh 



(18) 

The dimensionless frictional force F at the wall across 
one wavelength in the inclined channel is given by 



The pressure gradient is obtained from equation (14) 



F = Jh - 



0 



dp 

dx 



dx 



J 



dp _ (q+h)(2-N) 
dx (1-N) 



1 



„ oodinh Nsirhnh h „ t 
S, -1— +SJi 

m m m 3 



+r|ana 



f 



1 



_ cnhrrh Nsirhnh h ^ 
Sj -^+S>h 

V m m m 3 j 



ck+rpina 



(15) 



(19) 



The time averaged flow rate is 
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mi. Discussion Of The Results 

From equation (17), we have calculated the 
pressure difference as a function of Q for different 
values of coupling number N with a = 0.4, m=2, 

71 71 

rj = 2 and for different a (i.e a = 0, — ,— ) and is 

shown in figure (2-4). In figure( 3) It is observed that for 
chosen parameters the pumping curves intersect at a 
point in the first quadrant closure to q«o.18. This 

may be due to the inclination of the channel is 
horizontal. For Q<0.18 we observed that the 
pressure rise increases with the coupling number N. 

The behavior is otherwise when Q > 0.18. For free 
pumping the Q decreases with the increasing N. For 
free pumping we observe that the flux Q increases 



with increasing 



r 



a 



v 



0 < a < — 

2y 



The same behavior 



is observed for the pressure difference Ap when 
Q = 0. 



The variation of pressure rise with time 
averaged flow rate is calculated from equation (17) for 
different values of micropolar parameter 'm', and is 
shown in figures (5) to (7) for fixed a = 0.4, 

n = 0.2, T| = 2 .It is observed that the pumping 

curves that the pumping curves meet at a point 
between q = o.l and q = o.2 . this value is 

estimated as Q = 0.12, when q <o.i2 the 

pressure rise decreases with increasing m. The 
opposite behavior is noticed forQ > o.l 2 ■ For free 

pumping we observe that the flux Q increases with 
increasing inclination of the channel 



71 



0 < a < — 

2 j 



From equation (17), we have calculated 

pressure rise as a function of Q for different values of 

and is shown in figure (8) for fixed 



71 



\ 



a I 0 < a < — 

2 j 



that for an inclined channel, the peristaltic wave passing 
over the channel wall pumps against more pressure rise 
compared to a horizontal channel (a = 0)- For free 

pumping the flux Q increases with increasing angle of 
inclination a between 0 and^L ■ 

2 

The variation of frictional force with time 
averaged flow rate is calculated from equation (19) for 
different values of N, m, a for a fixed a = 0.4 and is 
shown in figure (9) to (13) it is observed that the 
frictional force F has the opposite behavior compared 
to pressure rise. For horizontal channel it is observed 
that the given N, the frictional force increases with the 

flux Q. 



Ap 



-1 

-2 
■3 



-4 — 
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Figure 2.The variation of Ap with Q for different values 

of n with a = 0.4, 171=2, r| = 2 , a =0 
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n - 0.2, T| — 2 , m - 2 and a = 0.4, ■ It is Observed Figure 3. The variation of Apwith Q for different values 
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of N with a = 0.4, m=2, T| = 2 , a = — 
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Figure 6.The variation of Ap with Q for different values 



Figure 4. The variation of Ap with Q for different values 



of m with a = 0.4, n=0.2, T| = 2 , a =Pi/4 
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of n with a = 0.4, m = 2 ,ri = 2 , cc=— 

2 





Figure 7.The variation of Ap with Q for different values 
of m with a = 0.4, n=0.2, T| = 2 , a = Pi/2 



Figure 5. The variation of Ap with Q for different values 



of m with a = 0.4, n=0.2, T| = 2 , OC=0 
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Figure 8. The variation of Ap with Q for different values 
of a with a = 0.4, 171=2, T| = 2 , n=0.2 



Figure 1 0.The variation of F with Q for different values 
of N with a = 0.4, 171=2, T| = 2 ,0C=Pi/2 





Figure 9. The variation of F with Q for different values 



Figurel 1 The variation of F with Q for different values 
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of N with a = 0.4, m=2, T| = 2 , a = — 

4 



of N with a = 0.4, 171=2, T| = 2 , OC=0 
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Figure 12. The variation of Ap with Q for different values 



of a with a = 0.4, 171=2, T| = 2 , n=0.2 
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Figure 13. The variation of F with Q for different values 

of m with a = 0.4, n=0.2, T| = 2 , 0C=Pi/2 
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